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Introduction 

A network is a graph with directed edges, with multiple copies of the edges 
allowed, and where each node v has an integer value specified, the so called 
excess of v, and each arc has an assigned positive integer value called its 
capacity. A feasible flow is an assignment of real values to the arcs of the 
network so that for any node v the difference between the sum of values in 
outgoing arcs minus the sum of values in incoming arcs equals the prescribed 
excess of the node v and the capacities of the arcs are not surpassed. In 
this paper we study the problem of effectively counting the number of differ- 
ent integral feasible flows in a network. It is well-known that this problem 
is #P-hard in the computational category of counting problems ^3] be- 
cause the problem of counting perfect matchings in bipartite graphs reduces 
to it. Despite this bad complexity, concrete applications abound in graph 
theory Jl], representation theory ^j], and statistics J2j and thus finding 
good methods for attacking concrete examples is of importance. Our goal 
is to show that using the algebraic-analytic structure of the problem allows 
us to count flows in complicated instances very fast, surpassing traditional 
exhaustive enumeration. Continuing the work started in [2] we present effec- 
tive counting algorithms from which one can in fact derive counting formulas 
when the excess function has parameters. 

The set of all feasible flows with given excess vector b and capacity vector 
c is a convex polytope, the well-known flow polytope, which is defined by the 
constraints §g x = b, < x < c, where &g denotes the node-arc incidence 
matrix of G (a network matrix). The incidence matrix $g has one column 
per arc and one row per node. Each column of &g has as many entries as 
nodes. For an arc going from i to j, its corresponding column has zeros 
everywhere except at the i-th and j-th. entries. The j-th entry, the head of 
the arrow, receives a —1 and the i-th entry, tail of the arrow, a 1. A famous 
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instance is the max- flow min-cut problem ^S]. This is the case when b has 
first entry v, last entry — v and elsewhere. In part (B) of Figure ^ we 
list all possible flows with v = 11, the maximal possible from the network 
information specified in part (A). 

For us, an important feature of the network incidence matrix &q is 
that it is unimodular. We say that the system $c is unimodular, if the 
columns of &g span a lattice, denoted by Z<3?g and, whenever a is in this 
lattice Z<3?g, the polytope P(&G, a ) = {x\ x > : <&g% = o} has vertices 
with integral coordinates. Even more strongly, network matrices are in fact 
totally unimodular matrices |19j . which means that the lattice generated 
by their columns is the standard integral lattice Z n . Note that the integral 
feasible flows are precisely the integer lattice points inside the flow polytope. 

Here is an example: The node-arc incidence matrix for the graph G\ in 
Figure El is defined by: 



The equation &Gi% = b reads as the series of equations x\ — x% = b\, 
X2 — %3 + x$ = 62, — x\ + X3 — X4 = 63. These 3 equations express the 
fact that, at each node v G {1,2,3}, the difference between the sum of 
values in outgoing arcs minus the sum of values in incoming arcs equals the 
prescribed excess bi of the node v. Feasible flows are restricted furthermore 
by the conditions < X{ < C{. 

The algorithm and formulas deduced here are based on the notion of total 
residue (see Section^), the main concept involved being the study of rational 
functions with poles on an arrangement of hyperplanes. The enumeration 
theory we present was extended to arbitrary rational polyhedra in [2S] . The 
particular description we do here is valid for all unimodular matrices (again, 
remember that a matrix A is unimodular if A has integral coefficients and 
the polytope ? = {if! M™|j4x = b}, has only integral vertices whenever b 
is in the lattice spanned by the columns of ^4). 

The following lemma implies that it is enough to describe our counting 
formulas and techniques for networks without restricted capacities on the 
arcs and that have no directed cycles; these are called acyclic uncapacitated 
networks: 

Lemma 1 Given a network G with n nodes and m arcs, with capacity c and 
excess function b, there is an acyclic uncapacitated network G with n + m 
nodes, 2m arcs, and excess function b (a linear combination of b,c) such 
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Figure 2: Network Gi with nodes 1, 2, 3, edges aci, X2, X3, £4, excess function 
b\ = 3,62 = —2,63 = —1 and capacity function c Xl = l,c X2 = l,c X3 = 
2; c Xi — !■ 



that the integral flows in both networks are in bijection. The network G is 
obtained from G by replacing each arc by two new arcs as illustrated in the 
figure below. 




Proof: For the network G with capacity c, the flows are the solutions of 
&G X = b, < x < c (*). There is a clear bijection (a projection) between 
the solutions of system (*) and the solutions of 
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The new enlarged matrix is denoted $g and called the extended network 
matrix. To the network G, with its set of nodes V and its set of arcs E, we 
have associated the new network G . The set of nodes of G is the disjoint 
union of the two sets V and E and the network G is obtained from G by 
replacing each arc by two new arcs as illustrated in the figure above: that 
is to each / 6 E is associated f\ = [/, j] and _/*2 = [/, i] where i is the tail 
of / and j is the head of /. Both arrows f\ and fi are oriented with their 
common tail {/} belonging to the set E and their heads {i} and {j} in the 
set V. Thus G is a directed graph, with n + m nodes and 2m arcs. If b G W 1 
is the excess vector and c € W 71 is the capacity vector of the network G, 
we define a new excess vector b G R n M. m . The projection of b on M n has 
coordinates i>, = 6j — YlfeE\taii(f)=i ca P ac ^y(f)- The projection of 6 of M m 
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is the capacity vector c. Let Tq be the matrix with one column per arc 
and one row per node defined as follows. The column corresponding to an 
arc has just one non zero entry: the tail of the arrow receives a 1. Then 
$G — Tq is the matrix with one column per arc and just the head of the 
arrow receives a —1. All other entries are 0. Thus 
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is equal by definition to the matrix the m first columns corresponding 
to new arrows f\, and the last columns corresponding to new arrows fa. 
Solutions of 
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are solutions of the equation 
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Thus we obtain a bijection between feasible flows of the network G with 
feasible flows of the uncapacitated network G. The correspondence assigns 
to the arc f\ the value xj and to the arc fa the value yf = Cf — Xf. 

Example 2 Consider the network G\ of Figure® Using the transformation 
of the previous lemma we would pass from the capacitated network to the 
uncapacitated network G2 illustrated in Figured and the excesses of its nodes 
are in the caption. 

Because of Lemma Q and due to interesting applications in represen- 
tation theory, it makes sense to focus our efforts on the special case of 
uncapacitated acyclic graphs, and we do so on Section |2j A particular case 
is what representation theorists would call the Kostant partition functions 
associated to the complete graph K n with n nodes. There are many ways 
to induce an acyclic orientation to the complete graph, here we take the 
following convention of orientation: whenever there is an edge of the graph 
G between i and j, with i < j, then we direct the arrow from i to j. 

Example 3 Consider the complete graph G on vertices 1,2,3,4. In this 
case, each vertex is joined to all the others and the incidence matrix of the 
network is 
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Figure 3: Network G2 with excess {2,-4,-3,1,1,2,1} and no capacities 
resulting from capacitated network in Figure |2J 



Another example of flow polytope is the Pitman-Stanley polytope |17j 
that is constructed starting from a multiple edge graph: 

Example 4 Consider the graph with vertices (1, . . . , n) and edges from {i, i+ 
1} and {i,n} and the last edge {n — l,n} of multiplicity two. In the case 
n = 3 then 



Also, within the class of flow polytopes, we will be investigating the 
famous transportation polytopes These polytopes are usually described 
in terms of m by n real matrices (denoted here by M TOn (R)): Fix c = 
(ci, . . . , c n ) e and d = (di, . . . , d m ) G R™ such that YT=i d i = Y^i=i c « 
and define T mtn (d,c) as the set 
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Then T m ^ n {d,c) is a polytope called the transportation polytope associ- 
ated for the vectors d, c. We can easily see that this is another flow polytope 
over a complete bipartite network K mn (see Figure [I] where the first m 
nodes receive excess values (d\, . . . , d m ) and the n nodes in the second block 
receive the excess values (— c\, — C2, . . . , — c n ). The arcs are oriented from the 
first block to the second. In the family of transportation polytopes there 
is a distinguished member, the Birkhoff polytope that has been extensively 
studied (see for instance the references in the recent paper (31). 




Figure 4: The transportation polytopes are network polytopes of complete 
bipartite graphs 

It is well-known that the counting formulas of integer flows in a network 
come in piecewise polynomial functions (see jSJIH]). It is therefore of interest 
to understand the regions of validity of each polynomial formula, the so 
called chambers. We dedicate in Section [3] some effort to understand the 
structure of the chambers and how to determine the number of chambers. 
The question of how many chambers are possible was first raised in |15j . 
The combinatorial investigations of the chambers for the partition functions 
was initiated by [Ij. See also [TUj . 

1 Formulas for the volume and the number of in- 
tegral points of flow polytopes. 

In this section, we outline the principles used in the algorithms we imple- 
mented for counting integer flows. The method is valid for general convex 
polytopes [21123], thus we describe things in a general setting when possible. 
In Section [21 we will use particular properties of flow polytopes associated 
with graphs to calculate the counting formulas. 

Let be an integral r by iV matrix with columns vectors <pi, ...,(/> n- Let 
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b be an r-dimensional column vector and V = {x e R+ |<I>x = b}, the rational 
convex polytope associated to $ and b. We assume that b is in the cone C(&) 
spanned by the non- negative linear combinations of columns (/>i,(/>2, ■ ■ ■ ,<f>N 
of $. Without loss of generality we may assume that rank($) = r. If this 
is not the case, take the subspace of W generated by the columns of our 
matrix and rewrite the polytope in term of an appropriate rank k matrix of 
dimension k by N. For example, for the network polytopes the matrices are 
not of full rank but deleting one of the rows turns them into one. 

In what follows we assume that kernel^) flR+ = {0}. Then is not in 
the convex hull of the vectors and the cone C(<J>) is an acute cone in R r . 
For a £ W we denote by 

P($,a) = {( Xl ,x 2 ,...,x N )€MN\ Zf=i^j = a}. 

It is obvious that P(<J>, o) is a convex polytope determined by the matrix 
<5. Define 

v($>,a) = volume(P(&, a)). 
If $ spans a lattice in R r and a belongs to this lattice, then define 

jfe($,o) = \P($,a) nz N \. 

Thus k($, a) is the number of solutions (x±, X2, ■ ■ ■ , xn), in non-negative 
integers Xj, of the equation J2f=i x j ( t ) j = a - The function a) is called the 
vector partition function associated to $. The name partition comes from 
the fact that if 3> = [ei, e±, . . . , ei] is the sequence of N times the standard 
basis vector of R, then P$(aei) n Z N is the set of solutions of the equation 
a = x\ + X2 + • • • + xn, that is the partition of the integer a in N integers. In 
particular, the function a — > k(&,a) depends strongly of the multiplicities 
in the system <3?. The basic starting observation is 

Theorem 5 Let z£l r denote a vector in the dual cone to C(3>). Then, 

1 



V k(^,a)e~^ = r-T. 



aeC($)nz r ii^e* 

The goal is to compute the inverses of these two equations. The point is 
that one can write efficient formulas for the inversion of Laplace transforms 
in terms of residues. In the sequel, we will write indifferently (cj), z) or 4>{z). 



8 



Let A + the set {«£}, this means the elements of $ are present without 
multiplicities. We define A = A + U — A + . A subset a of A + is called a 
basic subset if {a} form a vector space basis of M r . The chamber complex 
is the polyhedral subdivision of the cone C(A + ) which is defined as the 
common refinement of the simplicial cones C(o~) running over all possible 
basic subsets of A + . The pieces of this subdivision are called chambers. We 
will discuss the chambers in detail, specially how to compute the chambers, 
in Section |31 The important fact to remember is that for each chamber there 
is a quasipolynomial formula for a) and we explain now how to derive 
the formula on a given chamber. 

Each <j> E A determines a linear form on C r and a complex hyperplane 
{z € C r \(f)(z) = 0} in C r . Consider the hyperplane arrangement 

H C = |J{zEC r |0(*) = O} 

and let R/s. denote the space of rational functions of z G C r with poles 
on He- A function in i?A can be written -P(z)/n</>eA 4 > { z ) n4 ' where P is a 
polynomial function on r complex variables and are non negative integers. 
A subset a of A is called a basic subset of A, if the elements € a form 
a vector space basis for W . For such a, set 

After a linear change of coordinates, the function f a is simply 1 — 
and we denote by S<\ the subspace of i?A spanned by such "simple" elements 
f a . Elements f a are, in general, not linearly independent, as we see in the 
example below. 

Example 6 Let A + be the set A + = {ei,e2, (ei — C2)}- Then we have the 
linear relation 

1 _ 1 1 

xy y{x-y) x(x-y) 

between elements f^J^; f^z w ^ a i = { e ij e 2j-, 02 = {ei,(ei — e2)} and 
o~3 = {^2, (ei -62)} basic subsets of A + . 

Partial differentiation <9j preserves the space Ra- The key result we need 
is that there is a well-defined decomposition of R& under the action of partial 
differentiations, a free module part generated by the basic rational functions 
f a , and a torsion module part, which is unnecessary for calculations and can 
be neglected. 
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Theorem 7 (Brion-Vergne |9j) The vector space Sa is contained in the 
homogeneous component of degree —r of Ra and we have the direct sum 
decomposition 

r 
i=l 

We call the projection map 

TresA ■ Ra —> Sa 
according to this decomposition the total residue map. 

The projection Tres A (f) of a function / with poles on the union of hyper- 
planes Tic depends only of the smallest hyperplane arrangement 7i' c con- 
taining the poles of /. Therefore we just denote by Tres(f) the residue of 
a rational function / with denominator product of linear forms. 

Example 8 Observe that if we work in R 1 and A = {±ei}, then Ra is the 
space of Laurent series 

L = {/(*) = £ a k z k }. 

k>-q 

The total residue of a function f(z) £ L is the function —j-- The usual 
residue, denoted Res z= of, is the constant a_i- 

We denote by R A the obvious extension of R A , when we replace the 
space of polynomial functions on r variables by the space of formal power 
series on r variables. Let F : C r — > C r be an analytic map, such that 
-F(O) =0 and preserving each hyperplane <f> = 0. If / G Ra, the func- 
tion (F*f)(z) = f(F(z)) is again in Ra- Let Jac(F) be the Jacobian of 
the map F. The function Jac(F) is calculated as follows: write F(z) = 
(F 1 (z 1 ,Z2,---,z r ),...,F r (zi,z 2 ,...,z r )). Then Jac(F) = det((-£-Fj)ij). 
We assume Jac(F)(0) does not vanish. For any / in Ra, the following 
change of variable formula, which will be useful in our calculations later on, 
holds in Sa- 

Tres(f) = Tres(Jac(F)(F*f)). 

Note that the total residue of a rational function is again a rational func- 
tion. By definition, this function can be expressed as a linear combination 
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of the simple fractions f a (z). If / £ Sa, then Tres(f) is just equal to /. We 
also know that Tres vanishes on homogeneous rational functions of degree 
m, whenever m ^ — r and that Tres vanishes on derivatives. If / = |j ^ fc ^ 
(with P a polynomial in r variables) has a denominator product of linear 
forms (4>k, z ) which do not generate, then it is easy to see that / is a deriva- 
tive and the total residue of / is equal to 0. We are now ready to fix our 
notation and recall the key formulas. 

Definition 9 For a G R r , define 

e^ 1 (a,z) N ~ r 
,U(a)(z) =Tres(— N — -) = — ^Tres(- 



and its "periodic" version 

,(a,z) 



K$(a)(z) = Tres( 



nf =1 i- e -^) 



The equality: 



{a ' z) , 1 „ , (a,z) N - 



Tres( — tt ) = — — —Tres{ 



'Il£=i<&,*> (N-r)l ' V n£a<^.*> 
follows right away from the fact that the total residue vanishes on homoge- 
neous rational functions of degree m, whenever m ^ —r. 

By definition, J$(a)(z) and K$(a)(z) are rational functions of z homoge- 
neous in z of degree — r. They are polynomial functions of a of degree N — r 
and the homogeneous part in a of degree (N — r) in K${a)(z) is J$(a)(z). 

Example 10 Let us compute J$(a)(z) and K$(a){z) in the case of the 
Pitman- Stanley polytope associated to $g of Example ^| The matrix <3?g 
is a 3 by 4 matrix of rank 2 . Deleting the last row leads to 



110 
-10 11 



Then J$(ai, a 2 )(zi, z 2 ) = TVes ( ^i^^ 



is 



1 ^ ( {a\zi + a 2 Z2" 2 
-Tres 1 



2! V Oi - ^2)21^2 
-^-Tres{- 5- f ) + aia 2 Tres(- f ) + -^Tres{- ^ — 

2 (Zl - Z2)Z\Z% [Zl - Z2)Z\Z% 2 (Zi - Z 2 JZl 

— Tres{- — ^ ) + a\a 2 Tres\- — J H Tres{- 



\z\-z 2 )z\ ' {zi - z 2 )z 2 2 {z\-z 2 )z\ 
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Now 7 — ^ — and -, — ^ — are simple elements so that they are equal 

\z\—z 2 )z 2 (Z\—Z 2 )Z\ r a i 



Zl 



to their respective total residue. To compute the total residue of jz^T^xp > 
we write z\ as a linear combination of linear forms in the denominator, in 
order to reduce the degree of denominator: 

Zl (zi - z 2 )+ z 2 J_ 1 



(Zl - Z 2 )z\ (Zl - Z 2 )z\ z\ ' (Zl - Z 2 )z 2 

The total residue of \ is 0, as \ = is a derivative, thus Tresi-, — 21 , 2 

J 2j ' " z 2 22 y (Zl-Z 2 )Z 2 1 - 

t — . We finally obtain: 

\Z\—Z 2 )Z 2 J a 

\a\ + 2aia 2 1 a\ 
J<s>{ai,a 2 ){zi,z 2 ) = --, r V-- 



(1 - e-^- z '^){l - e- 2 i)(l - e" 22 ; 
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2 (zi - z 2 )z 2 2 (zi - z 2 ) zi ' 

We now compute: 

/ e (a 1 z 1 +a 2 z 2 ) 

K$(ai,a 2 )(zi,z 2 ) = Tres 
This is 

Tres ( 1 c (^i+a2z 2 ) (z 1 -z 2 ) z x z\ 

\(zi - z 2 )z x z\ (1 - e-^- 22 )) (1 - e~ z i) (1 - e" Z2 )V ' 

We replace the analytic function 

(a lZl +a 2 z 2 ) (fl ~ z l) £l z 2 

(1 _ e -(*i-»)) (i - e -zi) (1 _e-«) 2 

fry «is Taylor series at z\ = 0, z 2 = 0, and keep only its term N(ai,a 2 )(zi, z 2 ) 
of homogeneous degree 2 in z±, z 2 which is 

(t^ + ai + i ai 2 ) zi 2 + + a 2 + ^ ai + a\ a 2 ) z\ z 2 + a 2 + ^ a 2 2 ) z 2 2 . 

77ms i^$(ai, 02) (^i, Z2) is equal to 

Tres ( N(ai,a 2 )(zi,z 2 ) 

V (*i - ^2)^1^2 

Arguing as for J$, we finally obtain that K$(ai, a 2 )(zi, z 2 ) is equal to 
1 a 2 + 2aia2 + 3ai + 2a2 + 2^1 a 2 + a 2 



2 (zi - z 2 )^2 2 (zi - z 2 ) z\ 
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We are now ready to write the formulas to compute the volume and 
number of integral points. See [2J Section 2] for details. To each chamber c 
of the subdivision of C(A + ) is associated a linear form / — ► ((c, /}) on Sa- 
If the system $ is unimodular, as is the case for networks, it takes value 1 
or on fa whether or not c is contained in C(a). 

Theorem 11 (Baldoni-Vergne |2j) Let c be a chamber of the subdivision 
ofC(A+) 

1. For a£c, the volume of P(<&,a) is given by 

u($,a) = ((c, J*(a))). 

2. If the system $ is unimodular, then for a G c n Z3>, the number of 
integral points in -P(<J?, a) is given by 

k(<S>,a) = ((c,K*{a))). 

3. The function a \— * v(<&, a) is polynomial on a chamber c. 

4- If the system <3? is unimodular, as is the case for networks, the Ehrhart 
function a \— > a) is polynomial on a specified neighborhood of a 
chamber c. 

A more general formula for arbitrary $ spanning a lattice Z$ in W is 
given in [23] • Now, the question is how to apply these two formulas for 
the computations with flow polytopes. The calculation of total residues will 
simplify considerably. 

2 Counting Integer Flows in Networks 

In this section we will focus on flow polytopes for acyclically directed graphs. 
We already justified in the introduction this makes sense, as other networks 
can be reduced to acyclic uncapacitated networks. Consider a r + 1 real 
dimensional vector space. Let (the positive root system of A r ) be defined 
by 

At = {(e i -e j )\l<i<j <(r + l)}. 
Consider E r the vector space spanned by the elements (ej — ej), then 

E r = {a E M r+1 |a = a\e\+- ■ ■+a r e r +a T+ \e r+ i with ai+a2+- ■ -+a r +a r+ i = 0}. 



13 



The vector space E r is of dimension r and the map 



(1) 



/: M r 



defined by 



a = (01, 02, . . . , a r ) 1 — ► a = a\e\ + ■ • • + a r e r 



(01 + ■ • • + a r )e r+ i 



explicitly provides an isomorphism of with the Euclidean space W . Let, 
as before, <& = {(f>±, fa, • • • ,<Pn} denote a sequence of non-zero linear forms 
belonging to A$ . We assume that the vector space spanned by 3> is E r . 
This sequence is completely specified by the multiplicity my of the vector 
&i — ej in Explicitly for the transportation polytope T mjTl (d, c), if we 
denote by <J> m ,n C ^m+n-i the ro °ts associated to it, then we have <& m ,n = 
{(ej — e j ) 1 1 < i < m, m + 1 < j < m + n} and thus m-jj = 1 if 1 < i < 
m, m+l<j<m + n, m^j = otherwise. 

It is clear that the polytope P($,a) is the polytope associated to the 
uncapacitated network with (r + 1) nodes, where the arc i 1— > j (i < j) 
appears my times (my can be for some arcs), and with excess function Oj 
at each node 1, 2, . . . , r and — {0,1 + 0,2 + - ■ ■+a r ) at the last node r+1. Indeed 
we have seen in Remark 3 that the columns of the matrix corresponding to 
.P(3>, a) are vectors of the form — ej for some i and j. 

The hyperplane arrangement (setting z T+ \ = 0) generated by A + is given 
by the following set of hyperplanes: 



A function in R^ r is thus a rational function f(z\, Z2, ■ ■ ■ , z r ) on C r , with 
poles on the hyperplanes z; L = zj or Z{ = 0. The following result is proved 
by induction in |2j, Proposition 14. 

Lemma 12 Let S r be the set of permutations on {1, 2, . . . , r} and f n , f w ,w G 
E r 6e defined by 



{zi\\ < i < r} U {(^ - Zj-)|1 < i < j < r}. 



U{zx,z 2 ,...,z r 



) 



1 



{z\ - z 2 ){z 2 - Z 3 ) ■ ■ ■ (Zr-i - z r )z r 



and 



f w (z 1 , ...,z r ) = w-f w (zi, ...,z r ) 



1 



( z w(l) — z «)(2))( 2; w)(2) — z -. 



"m 



(3)) ' " " ( z w(r~l) ~ z w(r)) z w(r) 



then 
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(2) dim Sa t = r\ 
and 

(3) {f w (z 1 , . . . ,z r ) = w f 7T (z 1 , . . . ,z r ), WG^r} 
is a basis for Sa t 

The cone C(A^) generated by positive roots is the cone a\ > 0, 01 + 02 > 

0,. . . , a% + Gt2H \-a T > 0. We denote by c + the open set of C(A^) defined 

by 

c+ = {a G C(A r + ) | en > 0, i = 1, . . . , r}. 

It is a chamber of our subdivision, and will be called the nice chamber. 
The importance of this chamber is that its "permutations" form a "basis" 
for the formulas that express volume and number of integral points. If c 
is a chamber for C(<3?) then there exists a unique chamber of C(A+) that 
contains c. 

Definition 13 (|2j) Letrriij (i < j) be the multiplicity of the vector ei — ej 
in (i.e. this is the number of times the arc i,j is present in the network). 
Let N = £\ • n%i t j the total number of arcs. We explicitly write down the 
functions J$(a) and K§{a) for our choice of <£, a. Recalling that z r+ \ = 0, 
we have that 

• T*(n}(n 7 \ - 1 Trr-i I (a 1 z 1 +-+a r z r ) N - r \ 

• J$(U)y^l, . . . , 4, r ) — / TV — rV \ m l,r + l m 2,r + l m r,r+l j-, , xm^ )l 

V '' \ z l z 2 "' Z r lLl<i<j<r\ Z i~ Z j ) / 

• K *W*> ■ ■ ■ ' Zr) = Tres ( nL^^-^^C^a-^-^-^)^ ) • 

We now write these functions in two specific examples. 

Example 14 We consider the polytope associated to a complete bipartite 
graph with 3 nodes on each side. Recall that in this case the matrix that 
determines the polytope is given by the vectors $ = {e\ — e^e\ — e^,e\ — 
e&, e 2 - e 4 , e 2 - e 5 , e 2 - e 6 , e 3 - e 4 , e 3 - e 5 , e 3 - e 6 }. So 

( rriij = 1 1 < i < 3 and 4 < j < 6 

1 otherwise 

and 
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J$(a)(zi, . . . ,z 5 ) = \Tres 



(ai zi +Q2 22 +Q3 23 +"4 Z4+g5 £5 ) 4 

2l2 2 2 3 11 l<i<3 (^i-Zj) 
4<j<5 



i<r$(o)(zi,. . . ,z 5 ) = Tres 



e a l 2 l e a 2 z 2 e a 3 z 3 e a 4 z 4 e a 5 z 5 

[ l<i<3 1 
4<j<5 



Example 15 We consider the polytope determined by the complete graph 
K5, in other words = A\. We obtain 

. J* (a) (21, . . . , z 4 ) = ±Tres ( ("i^+^+^+^M 

TS f \f \ rr I e a l z l e tI 2 z 2e a 3 z 3e a 4 z 4 

In handling the formulas that we have for computing the volume and 
the number of integral points, the first problem is that of computing the 
total residue. This is in general a very difficult task. On the other hand, 
as we have seen, there is a very nice basis in Sa,. and this will allow us to 
rewrite the formulas in terms of iterated residue, which are certainly more 
tractable. The point is that one needs to find some, but not all, simplicial 
cones that contain the chamber determined by a. This is a step that allows 
the complexity of the algorithm to be reduce. We are now going to introduce 
the iterated residue for A r . 

Recall that, via the identification (^Q) of E T with W , a function in Ra t 
is a rational function f{z\, z 2 , ■ ■ ■ , z r ) on C, with poles on the hyperplanes 
Zi = Zj or Zi = 0. For a permutation a £ S r define the linear form on Ra v 

Ires° =0 f = Res z<1)=0 Res Zi7{2)= o ■ ■ ■ Res Za{r)=0 f{z 1 ,z 2 , ...,z r ) = 

Res Zl=0 Res Z2= o ■ ■ ■ Res Zr= of{z a -i^, z a -i^, • • • , z a -i^). 
In particular for a = id the linear form / \—> Ires z= of defined by 

Ires z=0 f 

= Res Zl=0 Res Z2=0 ■ ■ ■ Res Zr=0 f(z 1 , z 2 , . . . , z r ) 

is called the iterated residue. 
Remark 

• the linear form / Ires^ =0 f on Ra t induces a linear form on 5,4,,, 
since it vanishes on the vector space of derivatives Ya=i 9iRA r ■ 
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• Ires a z=Q f w = 5°. 

• the r! linear forms Ires° =0 f, a E S r , on Sa t are dual to the basis f w . 

Iterated residues are easier to understand, and we will see shortly how 
to use them in connection to our formulas. Let w € S r and n(w) be the 
number of elements i such that w(i) > w(i + 1) (this is called the number 
of descents of the permutation w in HO]). We denote by C+ C C(A r + ) the 
simplicial cone generated by the vectors 

e (l)( e «>(i)- e u>(2))) e(2)(e^(2)-e w ( 3 )), . . . , e(V-l)(e t( ,( r ._i)-e u ,( r .)), (e„,( r )-e r+ i), 

where e(i) is 1 or -1 depending whether w(i) < w(i + l) or not. When w = 1, 
then C\ = C(A^). The following lemma is easy to see. 

Lemma 16 Let a = ]Cj=i a j e j * n Er- The cone C E r is given by the 
following system of inequalities Y?j=i a w{j) — 0> f or a ^ * such that w{i) < 
w(i + 1), but Yl)=i a w(j) < if w(i) > w(i + 1). 

From Theorem we obtain: 

Theorem 17 (|2j) Let c be a chamber of C(<1>). Consider the set of ele- 
ments w £ S r such that c C C+ . Then, for f £ S 1 ^ , 

((c,/)}= £ (-i)^)/^^- 1 /. 

In particular for f = J$(a) we obtain 
Formula 1: /or a £ c, we /love 

«(*,a) = «c,J*(a)» = £ (-l)"H /res - =0 j $(a ). 

The formula is a direct consequence of the fact that Iresf =0 is the dual 
basis of f w . We have seen that to compute the number of integral points of 
our polytope we need to compute K$(a). Let tj = rrijj+i + • • • + rrij^+i — 1, 
where we recall that rriij is the multiplicity of the root ej — ej in <£. After a 
change of variable for the total residue, we obtain: 
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Theorem 18 Let a = ^1=1 a^i in E r n Z r+1 . Let 

(1 + zi) ai+tl (1 + z 2 ) a2+t2 •••(!+ z r ) a - +u 



f<$>(a)(z) mi, r+ i m 2 ,r+2 "V.r+i n . _. .,„,,• 

^2 ■••-Sr lll<i<j<A z « 

T/ien Formula 2: for a € c, 

fc(*,a)= ^ (-i) n(tu) ^r=o/*(«)- 

We now want to give an even more explicit formulation of the above result 
suited to be directly implemented. For this purpose we need to introduce 
some more notations. For a G E r , let def(a) be defined by def(a) = a + 
e SaG$ a + e2 (YJi=i e i - re r+i) with e = ^zr and m the maximum of the 
multiplicities m«. 

A wall of j4+ is a hyperplane generated by r — 1 linearly independent 
elements of A£ . The cells in C(A+) \ TL (TL being the set of hyperplanes 
for ) are open cells, interior of polyhedral cones. We will call these open 
cells topes. We will say that a € C(A^) is regular if a is not on any wall for 

. The walls of are easily characterized since they are the kernel of a 
linear form as ^igj a i where J is a subset of {1,2,..., r}. It is then easy to 
decide whether a vector a is regular or not. 

If a is a regular element we let c denote the unique chamber of C(A^) 
containing it. Then the set Sp(a) = {w £ E r |c C C+} can be computed 
without explicit knowledge of the chamber. In fact one can easily see that 
the set Sp(a) consists of those w £ S r that satisfy the following conditions: 

' if a w ( 1 ) > then w(l) < w{2) else w(l) > w{2) 

if a w m + a w (2) >0 then w{2) < u>(3) else w(2) > w(S) 

if a w(i) + • • • + <hu(i) — then w(i) < w(i + 1) else w(i) > w(i + 1) 

if a w n) + • • • + a u ,( r _ 1 ) > then w(r — 1) < w(r) else w{r — 1) > w(r) 

An element of Sp(a) will be called a special permutation. 

Remark that if cij > for all i < r, then a = Y^=i a i e i ~ (Ya=i a i) e r+i 
belongs to the closure of the nice chamber c + and Sp(a) = {id}. 

Now we can state Theorem UHl as follows: 



Theorem 19 Let & C be a system generating E r . Let a = Yli=i a « e « ^ 
E r , a r+ \ = — (ai + • • • + a r ), aj G Z and assume that a € C(A^). 
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Write 



(1 + zi) ai+tl (1 + z 2 ) a2+t2 ■■■(! + 2 r ) a ''+*'- 



/$(ai,a2, . . . ,a r )(z) — mir+ m2r . +2 m r ,r+i n , _ ^ 

z l z 2 ■ ■ ■ Z r lll<i<j<r\ z i- z j) 

Then 

• Formula 2A: if a is regular then 



k(*,a)= Yl i-l) n{w) Ires w z= ,U(a). 

■w€iSp(a) 

• Formula 2B: if a is not regular then 

fc(<M) = (-l) n{w) Ires™ =0 U(a). 

weSp(def(a)) 

Remark Formula 2B in the theorem follows by observing that the cham- 
ber containing the regular element def(a) contains a in its closure. The 
deformation has to be done with care to deal with some border cases. The 
following lemma, that we state for completeness, shows that the deformation 
with a,i integers is small enough to take care of such cases. 

Lemma 20 Given a <G C(A+) PI Z r+1 , define def(a) := a + eJ2ae® a + 

l 

2mr' 2 



e2 (S[=i e i ~ re r+i), e = 2mr' 2 w ^ ere m * s ^ e maximum of the multiplicities 
rriij . Then the following holds: 

• def(a) is regular, i.e. it belongs to a chamber. 

• if t is a tope and a £ r then def(a) £ r 

• a E C(A+) if and only ifdef(a) € C(A+) 

• In general if & is a subset of , a £ C(&) if and only if def(a) £ 
C(<&). 

For example, we obtain the following formula for the complete network 
K r+ \ on r + 1 nodes, with excess vector a\, a 2 , ■ ■ ■ , a r , a r+ \ = — Yli=i a i- m 
this case, the function k(A*~,a) is the so-called Kostant partition function 
and has special importance for the representation theory of the group GL(r+ 
1,C). 
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Corollary 21 For a G C(A^) n 27 +1 , i/ie Kostant partition function is 
given by: 

k(At,a )= V (-D-W/n^fi + Sg^Hll *>*7' -('+ *£ 



w/iere 

a if a is regular 
def(a) otherwise 

In particular, if > for 1 < i < r, we have 



i 

a 



Res Zl= oRes Z2= o ■ ■ ■ Res 



Z r =0 



k(A+,a) = 
[1 + zi) ai+r - 1 (l + z 2 ) a2+r - 2 •••(! + z r ) ar 



Z\- ■ ■ Z r Yll<i<j<r( Z i Z j) J 

Similarly we may write a formula for the transportation polytope T m ^ n (d, c). 



Corollary 22 Let a = Sj=i c j e m+j, with ali andcj non negative 

integers. Then the number of integral points in T m ^ n (d,c) is equal to 

£ (-ir^wf =0 

^ (1 + + Z2 )^+n-l . . . (1 + Zm )d m +n-l (1 + Zm+1 )-^i-l . . . (1 + Zm+w _ 1 )-'=n-i-l 

ZX ' • • Z m Y\ l<i<m (Z{ Z m ^-j) 
l<j<n-l 

, _ J a if a is regular 
|de/(a) otherwise 

2.1 The Algorithm for Counting Integral Flows. 

Scope of this section is a brief description of the various algorithmic proce- 
dures that were implemented with the symbolic language Maple and that 
achieve the formula for the number of integral points described in Theo- 
rem^! This software is available at www. math. ucdavis . edu/~totalresidue. 
The initial data are an r by N matrix A whose columns are the elements of 
<1> and an element a = {ax, ■ ■ ■ , a r } S Z r that determines the polytope. The 
ingredients that we need to compute are: 
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1. The element a' = def(a) obtained by deforming the initial parameter 
a. 

2. The set of permutations that appear in the formula, that is the set of 
special permutations Sp(a'). 

3. The residues that appear in Formula 2. 

We will discuss the ingredients for each one of these steps listing the 
various algorithms that are related to the part we are describing. 

First of all we want to check if our vector is in C(A£), that is in the cone 
generated by {(ei — e<i), {&i — 63), . . . , [e.( r -\\ — £r),&r} because otherwise the 
polytope is empty and there is nothing to do. To be in the cone, a must 
satisfy a\ > 0, a%+a2 > 0, . . . , 01 + 02 + - ■ - + a r > 0. The procedure check- 
vector verifies whether this is true or not. In fact because of Lemma 1201 
we may use def(a) instead of o and we do this to simplify the procedures. 
We compute the element def(a) via the Maple procedure def-vector. The 
vector def(a) is used in all the formulas defining Sp(a) instead of a, whether 
or not a is regular. This takes care of the first part. 

For finding the subset Sp(a) of E r , we use the procedure special-permutations. 
We stress that using the Maple function combinat[permute] is impracti- 
cal and does not go very far because of memory limitations. Our ap- 
proach constructs recursively the permutations subject to our conditions, 
thus we save much memory in listing only those permutations. The set 
Sp(a) depends strongly on the element a. We do not have upper bound 
estimates on the subset Sp(a) C E r; but it seems that this set is small 
compared to E r . One of the worst experimental cases for the complete 
graph Kio on 10 nodes (the case of Ag) is the case of the vector a = 
[30201, 59791, 70017, 41731, 58270, -81016, -68993, -47000, -43001, -20000] 
where the number Sp(a) C Eg is 9572, certainly much smaller that 9!. Ex- 
periments show that the time spent to compute this set is rather small. 

Each permutation it; £ Sp(a) gives rise to the simplicial cone con- 
taining a, this corresponds to a vertex of the polytope P(A+,a). However, 
clearly the cardinality of Sp(a) is much smaller that the number of vertices of 
the partition polytope P(A+,a). For example, for a = [a%, 0,2, ■ ■ ■ , a r , — QZi=i Oi)] 
with dj > 0, we have already remarked that the cardinality of Sp(a) is 1, as 
Sp(a) is reduced to the identity permutation. 

Finally, for the last step we need to compute the residue. Recall that we 
need to compute 
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w (1 + Zi) ai+tl (l + Z 2 ) a2+t2 •••(! + ZrY^r 
2=0 mi,r+l m 2 , r+2 mr, r +l TT f r _ , ]">,, 

z l Z 2 ■■■£■{■ lll<i<i<rV^* z 3> 

with w one of the special permutations. Let us denote by F the function 
appearing in the formula above. The function F is a product of a certain 
number of functions. This allows us to take the residues by introducing 
little by little the part of the function F containing the needed variable. To 
make things simpler we assume that w is the identity permutation. We start 

by taking the residue at z r = of the function q := mr r li „Z 1 , -z— 

Z r T ' r+ WjZlizj-Zr) "V 

Suppose g r (z±, Z2, ■ ■ ■ , z r -±) is the result. We continue by taking the residue 
in z r _i of the function g r multiplied by all the factors of the original func- 
tion F that involve the variables z r _i and so on. The way we compute 
the residue in one variable z of a function g(z) = F(z)/z u , where F is 
analytic, is by computing the Taylor expansion of F up to the estimate 
we have for the order u of the pole of the function g and then taking the 
coefficient of 1/z. The argument just described is implemented via differ- 
ent procedures: coeex,invi,trunc-next-function and RRK . Finally, the 
procedure numb er-kost ant adds up, with a sign (the appropriate sign is 
computed using segnop), all residues coming from the different special per- 
mutations, thus getting Formula 2. The procedure polynomial-kostant 
computes the polynomial a t— > k(&, a) on the chamber determined by a. 

As we pointed out we need an uniform estimate for the order of poles 
appearing. The result for the order of pole is the content of the subsection 
that follows and it is implemented in procedure E. 



2.2 Estimates for the order of poles 

Let G r be a Laurent polynomial in the r variables z = (z±, z 2 , ■ ■ ■ , z r ) and 
let D r = r[i<i<j<r( z * — z j)- We have seen that we need to compute iterated 
residues of the form : 



Res Zl=0 Res Z2= o ■ ■ ■ Res Zr=0 G r /D™. 

The following key lemma will handle the situations that will appear in com- 
puting the estimate we are looking for. 

Lemma 23 Assume that G r = f^ zl '"'' Zr ? a H r (—, ...,—) where F is analytic 
and H r is a homogeneous polynomial of degree h , then 

Res Zr=0 G r /D™ 
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is a linear combination of functions of the form G r -\/ D™^ with 

F( Zl , . . . , z r -i) u 1 1 

W-l — 7 \(n+m\ U r-l{ , ■ ■ ■ J ) 

where H r -\ is a homogeneous polynomial of degree at most g + h — 1 and 
F(z±, . . . , Zr-i) is analytic. 

Proof: Let us prove the lemma for a monomial H r = z\ x ■ ■ ■ z^Z\Z^ where 
ii, 12, • • • , i r are non- negative integers such that i\ + i 2 + • ■ ■ + i r = h. We 

write Ill<i<r-l(^ " Z rT = (*1*2 • • • Zr-i)™ IWr-lt 1 ~ tJ 

The Taylor expansion of tt 1 n _z I1 w at z r = is 

lll<z<r— 1\ z ■ / 

E -|C/l| -|t/ 2 | -\Ur-A _|i7l| + |{72|+-+|J7 r _l| 

z 1 z 2 z r _ 1 z r 

C/i,..,f/ r _i 

where J7 S = {jf , jf , . . . , j^} varies over the m tuples of non negative integers. 
Write also F(z±, . . . , z r ) = ^ fc . . . , z r _i)z^. Thus we obtain 

Res Zr=Q ^ = 
-Res Zr=0 - 



7~ bl ■ ■ ■ 7 V_1 i 

( Z l Z r-1 1 W 

y (z 1 Z 2 ---Zr-l) 9+m A"V 

£ (-Pfc(^i, . . . , z r -i) z -^\ . . . z ;\ u r A ) 

k=0 U 1 ,...,U r -i:\U 1 \+-+\U r -i\=g-l+i r -k 

For 0<k<i r +g — 1, the monomial 

_-<i -v-i -|t/i| -|C/ 2 | -|f/r-i| 

z l z r- 1 z l z 2 z r-l 

is such that 

ii -\ Mr- 1 + l^i H 1- | C/, — x | = ii H hv-i + ir + g- 1 -k < h + g — 1 

and we obtain the lemma. 

Observe that if F = 1 then the same proof shows that i7 r is homogeneous 
of degree precisely h + g — 1. Now starting from G r = ^f^jy^r we want to 
compute 
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Res Zk+1=0 Res Zk+2=0 ■ ■ ■ Res Zr _ 1=0 Res Zr=0 G r /D™. 
Applying the lemma with h = 0, we obtain that 

Res Zr=0 G r /D r r 



r 



G 

is a linear combination of functions of the form where 



F(zi,...,z r -i) , 1 1 
G r-i - -r ^—H{—,... 



[Z\ ■ ■ ■ Z r -i) 2m 01 ' " ' ' Zr-i 

and -ff is homogeneous of degree at most m — 1 , thus at the next residue we 
get again a linear combination of functions of the form n m 2 where 

_ F( Zl ,...,z r - 2 ) „. 1 1 v 

(Zl ■ ■ ■ Z r -2r m Z\ Zr-2 

with H homogeneous of degree at most 2m + m — 1 — 1 = 3m — 2, so finally 
the last residue in z^+i = leaves a linear combination of functions of the 
form 

Gk 

with 

F(zi,...,z fc ) ^ 1 J_ 
(zi---z fc )(«-*+i)»» ai zi'-'V" 
Here H is homogeneous of degree at most ( r fc )( r ^ fc+i)m _ ^ _ ^ j n p ar _ 
ticular, considering H(—,... we have the estimate on poles we were 
looking for. 

Corollary 24 1. Let G r = F fa.'.' z '\m , with F analytic. Then the function 
Res Zk+1=0 Res Zk+2=0 ■ ■ ■ Res Zr _ 1=0 Res Zr=0 G r /D™ 
has a pole in Zf. of order at most m ( r fc )( r _ ^ r _ ^ 



2. In particular with the notation as in Theorem \18l if 
then the pole in a{zk) of the function 



m = maximum* ,m 



Res z a <.k+i)=a ■ ■ ■ Res Zcr(r)=0 U(ai,a 2 , a r )(z) = 

(1 + Zl ) ai+tl (l + z 2 ) Q2+t2 ■■■(! + z r ) a - +t - 

z l z 2 ■ ■ ■ Z r lll<i<j<r\ z i Z 3I 

has at most order m ( r fc )( r _ ( r _ ^ independently from a G S, 
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3 The Chamber Complex 



In this section we discuss the chambers and how to compute them. It is 
important to emphasize that everything that we present in this section is 
valid for general matrices, not necessarily unimodular. There is an im- 
plementation of these ideas in the Maple program chambers available at 
www.math.ucdavis.edu/~totalresidue. Let A + the set of distinct vectors 
{<!>}. Recall the chamber complex is the polyhedral subdivision of the cone 
C(A + ) of nonnegative linear combinations of A + . It is defined as the com- 
mon refinement of the simplicial cones C(o~) running over all possible basic 
subsets a of A + . To be more precise we introduce now notation and the key 
definitions. In what follows, when we consider a subset I = {si, 82, ■ ■ ■ , Sfc}, 
where the elements of / are subsets of a set X, we assume there is a 
partial order on / by containment. Thus the set of minimal elements of I is 
denoted by minimalize(I) . We adopt the convention that the intersection 
of an empty family of subsets of X is X itself. 

Let A + be the set {(pi, 4>2, ■ ■ ■ , <Pn} of vectors in W . Recall that a wall is a 
hyperplane in lR r spanned by (r — 1) vectors of A + . Each wall W partitions 
the set of indices {1,2, ...,N} into three sets: zeros(W) = {i\4>i G W}, 
and two disjoint subsets posiW), neg(W) whose union posiW) U neg(W) 
is precisely the subset of {1,2, ... ,N} \ zerosiW). We consider the set 
{posiW), negiW)} = {neg(W) , pos(W)} . We denote by B the set of subsets 
a of {1,2, ...,N} such that a is of cardinality r and the set of vectors 
{4>i\i G a} are linearly independent. For convenience, we continue to call 
such a a basic subset of A + , thinking of a as a subset of integers or as a 
subset of elements of $ labeled by indices. 

For a G B, we consider the closed cone C(a) generated by a. If I is a 
subset of B, let F(I) = n ae iC(a) be the intersection of the cones C(<r), when 
a runs in /. We will say that / is a feasible subset of B if the interior of F(I) 
is non empty. A combinatorial chamber I is a maximal feasible subset of B. 
The polyhedral cone F(I) will be called a geometric chamber. The actual 
chamber Chamber(I) is the interior of F(I). Reciprocally, the collection 
/ is entirely determined by F(I). We have I = {a G B\F(I) C C(a)}. 
The collection of all geometric chambers and their faces forms a polyhedral 
complex that partitions the cone C(A + ), the so called chamber complex 

Unsung. 

Figure |S] shows an example, the chamber complex for the cone associated 
to the acyclic complete graph K4 we discussed in the previous section. The 
picture represents a 2-dimensional slice of the cone decomposition (the cone 
is 3-dimensional and pointed at the origin). The 6 dots labeled (e, — ej) on 
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Figure 5: A slice of the chamber complex for K4 



the drawing are the intersections of the rays M + (ej — ej) with the hyperplane 
(3xi+X2— X3— 3x4) = 2. Seven chambers, numbered from 1 to 7, are present. 
In the configuration of vectors of Figure 03 there are seven walls, one for each 
of the distinct lines obtained from the vectors in the configuration. 

Let 7i denote the hyperplane arrangement consisting of all walls. 7i 
contains as a subset the walls of the chambers. The cells in C(A + ) \ 7i 
are open cells, interior of polyhedral cones. We will call these open cells 
topes (following the oriented matroid terminology 7 ). Note that the set of 
topes is (typically) a much finer subdivision of C(A + ) than its chambers. 
See Figure El for a comparison between the chamber complex and the tope 
complex of the hyperplane arrangement Tt associated with the example in 
Figure 03 



e-e 3 e 2 -e 3 




e-e e , _e e-e 

1 2 4 1 2 



Figure 6: 8 topes (left) versus 7 chambers (right) 

A tope t of C(A + ) does not touch any wall of A + . Then, for each 
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wall W, we denote by posiW, r) the set of elements i € {1,2,..., N} such 
that fa £ A + lies on the same open half-space determined by W than the 
tope r. We say that posiW, r) is a non-face (this terminology is justified 
because these are the non-faces of a certain simplicial complex in the sense 
of Chapter two of [U). We denote by Chamber {t) the chamber containing 
the tope r. 

To each tope r, we associate the family of positive non-faces determined 
by the tope r (we have a non-face for each wall). Let us call this full 
family Polarized(r). Consider the family MNF{t) of minimal elements 
of Polar ized(r). This is the family MNF{t) = minimalize{Polarized{T)). 
The first main observation is that we can reconstruct the chamber Chamber (t) 
containing the tope r from the set MNF{t). This is very useful to construct 
one initial chamber. Later all others will be found from it. 

The set MNF(t) is a set of non-faces. Let / be the cardinality of the 
set MNF(t). Let us list MNF(t) := {pi,P2, ■ ■ ■ ,Pf}- Each pi is a non- 
face. We construct the family V(t) of sets v of the form v := {i\, 12, ■ ■ ■ ,if} 
with ii G pi, 12 £ P2-,- ■ ■ , if € Pf. These we call transversals of a family 
of sets. This family is denoted by transversal{MNF(r)) in the computer 
program we present. Again V(t) is a set whose elements are sets of indexes, 
its elements being subsets of {1, 2, ... , N}. The cardinality of a set v E V(j) 
may be smaller than /, as the family MNF(r) does not consists of disjoints 
sets. It is important to observe that if v is in V(t), then for any wall W, 
the intersection v (~)pos(W, r) is not empty. We have the theorem. 

Theorem 25 The minimal elements of the family V(t) := transversal(MNF(r)) 
are exactly the basic subsets a of A + such that r C C(o~). 

In other words, given the set MNF(t) associated to a tope r, the family 
of basic subsets a of A + such that r is contained in C{a) is precisely the set 
minimalize(transversal(MNF(T))). We are going to prove this theorem. 
We start by a lemma. 

Lemma 26 Every v € V{t) is such that the set of vectors {(f>%\i S v} gen- 
erates W . 

Proof: Let us see that a set v £ V{t) generates W . Indeed, if not, 
the set of vectors {_<j>%\i € v} would be contained in a wall W. Consider 
the set pos(W,r) and a minimal element p of the family MNF(t) : = 
minimalize{Polarized{r)) contained in pos(W,r). Then p (meaning the 
set of elements fa indexed by p) is contained in one of the open half-space 
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determined by W. Thus, contrary to our hypothesis, we would have vC\p = 0. 
QED 

We go on proving Theorem 1251 
Proof: Let a be a basic subset of A + (a (elements indexed by a) generates 
a simplicial cone). We now prove that if r C C(a), then a G V{t) and is a 
minimal element in the family of tranversal sets V(t). 

For each wall W, the set a Pipos(W, r) is non empty. Otherwise a would 
be contained in the closed half space determined by W, but would be on 
the opposite to r with respect to W, and the cone C(a) will not contain r. 
Let us pick for each p G MNF(r) an element <p p G a dp. It follows that a 
contains necessarily the set v := {4>p\4>p 6 °~ (lp;p G MNF(t)}, belonging 
to the family V{t). But then a = u, as a is a basic subset of A + and v 
indexes a set of generators of W by Lemma l26l Furthermore a is minimal, 
as all sets belonging to the family V(t) have cardinality at least equal to r. 

We now prove the converse. Let v be a minimal set of V(t). We claim 
that r is contained in the cone C{y). Otherwise, there would be a wall W 
separating r and C{v). But by construction of v there is an element 
contained in pos(W, r); a contradiction with W separating C(v) and r. Now 
all we have to prove is that v has cardinality r. 

Let x be a point in r. By Caratheodory theorem, there is a basic subset 
a contained in v such that x G C{a). Then the tope r is entirely contained 
in C(cr) because a tope is, by definition, not separated in two by any hyper- 
plane. The set a belongs to V(t) by the preceding discussion. But a C v 
and v is minimal, thus v = a. 

So we conclude that the set Chamber (r) of basic subsets er of A + such 
that r C C(cr) is the set minimalize(V(T)) of minimal elements of V(t) = 
transversal(MNF(T)). QED 

The lexicographic tope is the tope containing the vector £ = </>i + e02 + 
e 2 </>3 + • • • where e is a small number. The lexicographic chamber is the 
chamber that contains the lexicographic tope. 

Corollary 27 The following algorithm determines the r-simplicial cones 
C{a) that contain the lexicographic chamber associated with a particular 
labeling of the elements of A + , by finding the basic sets a that define them. 

1. Create the list L of lexicographic nonfaces pos(W,r) where r is the 
lexicographic tope, and W runs over all possible walls of A + . 

2. Let F = {A\, A2, ■ ■ ■ , A m } be the minimal non-faces from L. 

3. Find the transversal sets to the family F then minimalize the set of 
transversals. The result is o~i, . . . ,o~k the desired basic sets. 
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Now we are concerned with producing all other chambers from one initial 
chamber, such as the lexicographic chamber. For this we need to understand 
the polyhedron F(I). This is a pointed polyhedral cone. We recall, say from 
Chapter 8 in the book [T^J, that for a polyhedron P (e.g. F(I)) given by 
a finite set of inequalities Ax < b, a supporting hyperplane is an affme hy- 
perplane {x\cx = d} such that d = max{cx\Ax < b}. A subset of P is a 
face if F = P or F is the intersection of P with a supporting hyperplane of 
P. A facet of P is a maximal face distinct from P. We say a wall W is an 
essential wall of the geometric chamber F(I), if F(I) D W is a facet of the 
pointed polyhedral cone F(I). This is equivalent to W being a supporting 
hyperplane of F(I) and dim(F(I) D W) = r — 1. We say that two geometric 
chambers and F(I') are W -adjacent if they share a common essential 
wall W and dim(F(I) n -F(i"') H W) = r — 1. In particular, the wall is an 
interior wall. In what follows, unless is necessary to avoid ambiguity, we will 
simply refer to "adjacent chambers" without specifying the wall they share. 
We present now an operation that allows us to move, under certain con- 
ditions, from a geometric chamber to another adjacent geometric chamber. 
Since the geometric chambers form a connected polyhedral complex, we can 
then apply some standard search procedure, such as depth-first search, to 
enumerate and list all chambers. 

We denote by W the set of subsets v of {1,2,..., N} such that v is of 
cardinality r — 1 and the set of vectors {4>i\i G v} are linearly independent. 
In other words, if v is in W, the vector space C(v) spanned by the vectors 
{<j>i\i G v} is a wall W. If W is a wall we denote by W(W) the subset of W 
with elements those v such that C{v) = W. 

If v is in W, we consider the subsets zeros{v), pos{v) and neg{v). If i 
is not in zeros(u) , then v U {i} is an element of B. We denote by 
the subset of B consisting of elements a = v U {i} where i runs in pos{v); 
denote 5~{v) the subset of B consisting of elements a = vU{i} where i runs 
in neg(y)\ 

If W is a wall, and o a subset of {1,2,..., N} we denote by a n W = 
<7 n zeros(W) . We denote by B(W\facet) the subset of £? consisting of those 
elements cr such that a n is of cardinality (r — 1). In other words, W is 
spanned by a facet of the cone C{a). We denote by B(W\cut) the subset of 
5 consisting of elements a such that both sets a (lpos(W) and a D negiW) 
are non empty. For any subset I of we denote by 7(W|/acet) = / D 
B{W\facet) and by J(W[cut) = I D B(W|cut). 

Let 7 be a combinatorial chamber which is a maximal feasible subset of 
B . Let W be a wall, we define B(W,I) = {a n W[<7 G J(W[/acet)}. This 
is a subset of W(W) = {v G W|£(v) = VF}. If is an essential wall of 
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F(I), then (as we will see later) for each subset v 6 B(W, I) either 5 + {v) is 
contained in / or 8~(v) is contained in /, but not both. 





W 






plus 


\ \ minus 



Figure 7: A reflexion exchanges the simplicial cones supported on opposite 
sides of a wall. 

If W is an interior wall then define the reflexion operation, this is a new 
combinatorial chamber denoted by reflexion(I , W). We keep in reflexion(I, W) 
all elements a E I(W\cut), while we replace each subset 5 + (v) C I(W\facet) 
by its opposite 5~{v). The operation of reflexion has also received the name 
of flip by several authors. Applying a reflexion over any wall may not yield 
an adjacent chamber, as we see in the example of Figure El 



2 3 




6 5 



Figure 8: A reflexion using the wall 1,4 does not give a chamber 

The important fact is that if one performs the reflexions over essential 
walls the result is the desired one: 

Lemma 28 IfW is an essential interior wall ofF(I), and let reflexion(I , W) 
the geometric chamber obtained by reflexion of I along the essential wall W 
. Then the set reflexion(I, W) is the combinatorial chamber associated to 
the W -adjacent chamber sharing W with F(I). 
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Clearly all elements a € I(W\cut) and elements in 8~(u), when v runs 
over B(W, I), give rise to simplicial cones containing the W^-adjacent cham- 
ber. Conversely, any a in B such that the cone C(o~) contains the VF-adjacent 
chamber is either in I(W\cut) or in a set of the form 5~(v), with v £ B(W, I). 

The above lemma stresses the importance of determining the essential 
walls and that is what we describe next. Each essential wall W is described 
by a linear inequality, that reaches equality at F (I) P\W. The chamber is 
contained in the corresponding half-space. The presentation we have of the 
chamber is as the intersection of simplicial cones, their facets provide us 
with a system of inequalities whose solution is precisely the chamber. The 
trouble is that this system contains redundant inequalities. An inequality is 
redundant if it is implied by the other constraints in the system, so redundant 
inequalities can be removed. 

Our algorithm for finding the essential walls is based in the following 
statement, which is essentially Theorem 8.1 in page 101 of |19| . Here we state 
it for full-dimensional polyhedra (thus no equality constraints are present): 

Theorem 29 If no inequality in the system Ax < b defining the full- dimensional 
polyhedron P is redundant, then there exists a one-to-one correspondence 
between the facets of a polyhedron and the inequalities in Ax < b given by 
F = {x € P\aiX = (3i\, for any facet F of P and any inequality aix < (3i 
from the system Ax < b. 

So if we manage to remove redundant inequalities from the original sys- 
tem of inequalities associated to F(I) we would have found the essential 
facets of the pointed polyhedral cone F(I). To do this let us describe a 
direct method. Let Ax < b, s T x < t be a given system of m + 1-inequalities 
in d-variables x = (xi,X2, ■ ■ ■ , Xd) T ■ We want to test whether the subsystem 
of first m inequalities Ax < b implies the last inequality s T x < t. If so, the 
inequality s T x < t is redundant and can be removed from the system. A 
linear programming formulation of this is rather simple: 

/* = maximize s T x 

subject to Ax < b 

s T x < t + 1. 

Then the inequality s T x < t is redundant if and only if the optimal value 
/* is less than or equal to t. By successively solving this LP for each untested 
inequality against the remaining inequalities, one would finally obtain an 
equivalent non-redundant system. Thus the algorithm to recover all the 
essential walls as follows. 
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1. Find the inequalities of each of the simplicial cones in F(I). 

2. Remove redundant inequalities using linear programming until there 
is no redundant inequality left. By the previous theorem the wall is 
uniquely determined by setting to equality the inequalities. 

Thus to find all the chambers, we have 

Corollary 30 The following algorithm finds all the chambers of the vector 
set A+; 

1. Find the lexicographic chamber Lmitiai- Put that as the first element of 
a list of chambers L. 

2. Pick an element I of L for which we have not yet found its adjacent 
chambers. Determine its essential walls W using the method above. 

3. Perform the reflexions reflexion(I, W) = I(W) for each essential 
interior wall W . 

4- Add the I(W) to the list L of existing chambers if not already there, 
and continue until we have found adjacent chambers for all elements 
in L. 

Although we have a concrete algorithm now to generate all chambers for 
practical reasons it is highly desirable to improve the speed on recognizing 
the essential walls. For this we prove some necessary conditions of the 
essential walls of a chamber: 

Proposition 31 Let I be a combinatorial chamber (a maximal feasible sub- 
set of B). Let W be a wall of A + . If W is an essential wall of F(I), then 
the following conditions hold true: 

1. I = I(W\ facet) U I(W\cut). 

2. I{W\facet) + 0. 

3. For each v E W, either 

S+(u) n / + 0. Then 5+{v) C / and 5~{y) D I = 0; or 5~{v) n I + 0. 
Then b~(y) C / and 5 + {v) n I = 0; 

o 

4- Assume I(W\cut) is not empty . Then f\ a ei{W\cut) intersects W 

in an (r — 1) dimensional set. 
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We start the proof. Let {v±, . . . , v r -±} be independent vectors in W, 
generating a cone contained in F(I) (~l W. If a = {<f)±, . . . ,(f) r } £ I we 
denote by A a the matrix expressing {v±, . . . ,v r -±} in terms of a, that is 
v i = Sj=i a ji4>j- The matrix ay t has non negative entries for any a £ I. 

Denote by 

A,ct = [&i,l> • • • j &i,r— l] 

the components of vi,V2,-- - ,v r -\ on fa. These are the columns vectors of 

At- 

Lemma 32 Assume W is an essential wall of F(I). Suppose W is spanned 
by the vector set {v±, . . . , v r -±}. Then for each a £ I, either 

a) A^ c / for all i, 
or 

b) there exists an index k such that A^^ = while A Sy(J ^ 0,s ^ k. 

If a verifies the condition a), then a £ I(W\cut). If a verifies the con- 
dition b), then a € I(W\facet). 

Proof: 

Indeed, suppose that by rearranging the indices A\ a = 0, . . . , A qa = 
0, then the vectors {v±, . . . , v r -±} belong to the linear space spanned by 
{(j)q+i, . . . , fa} forcing r — q > r — 1 that is q = 1. Thus if a G a verifies 
either a) or b). 

Suppose we are in the first case. We now prove that a G I(W\cut). Let 
us see that if all the Ai )(T are non zero vectors , then there exists an element 

X eWn C(a) which would force that a G I(W\cut). Let X = YJlZl fm 
with fi > 0, then X G W. On the other hand X = fi Ej=i a ji4>3 = 

Sj=i(S!=i fi a ji)4>3 = Sj=i bj<frj- Because all the vectors A,ct are nonzero, 
then bj > for all j, thus X belongs to the open simplicial cone spanned by 
a. 

Consider the second case. Suppose for simplicity that iq = 1. Then 
{v\, . . . , v r -i} is a subset of the linear span £{(p2, ■ ■ ■ , fa}, therefore W = 
£{fa, ■ ■ ■ , fa} and hence a £ I(W\facet). 

Lemma 33 Let W be an essential wall ofF(I). Then the relative interior of 

o 

the cone generated by F(I)(~)W (inW ) is contained in f\ a ei{W\cut) 
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Proof: For X in the relative interior of the cone F(I) n W , we choose, as 
in the preceding proof, {v±, . . . , ty_i} independent vectors in W , generating 
a cone contained in F(I) n IF and such that X = Y7i=i fi v i> with fa > 0. 
Arguing as before, we see that X is in the interior of C(o~) for any a G 
J(W|cut). 

Lemma 34 Lei IF 6e an essential wall of F{I). Then I(W\ facet) is not 
empty. 

Proof: If I (W\ facet) was empty, then the open chamber Chamber(I) 

o 

would be Chamber (I) = ^ a &i(w\cut) C{a) and would intersect W due to 
the preceding lemma, and would not be contained on a half-space of W. 

If W is an essential wall of F(I), then F(I) is on one side of W, thus we 
have a distinguished non face pos(W,I). For each v G W spanning W, we 
have a distinguished set 5 + {v) = {v U {i}\ i G pos(W, I)} of elements of £>, 
while b~~{v) = \y U {i}| i G neg(W, I)} is disjoint from J: 

Lemma 35 If W is an essential wall of F(I), then for every v G W suc/t 
i/iai £(i/) = IF, t/ien 

1) We have <J-(i/)n/ = 

If5 + (v)nl ^ 0, taen C J. 

Proof: Condition 1 is clear, otherwise F(I) would be on the wrong side 
of IF. Now let x G Chamber(I) very closed to IF n F(I), and X in the 
interior of IF n F(I). Assume that a = v U {i} belongs to /. Then the point 
X is in the interior of C(u). The line [x, X] is in the chamber Chamber(I) 
except at the last point X. It cannot cross any boundary of any simplicial 
cone. Thus we see that it stays entirely in the interior of any simplicial cone 
spanned by v and a vector <j)j. with on the same side than x, as clearly 
its beginning (1 — e)X + ex is inside this simplicial cone. 

Thus we have proven that if IF is an essential wall of F(I), the wall IF 
satisfies 1, 2, 3, 4 in the statement of Proposition 1311 

Corollary 36 IfW is a wall of F(I) satisfying 1), 2), 3) and not 4) ; then 
reflexion(I ,W) is not a feasible subset ofB. 

Proof: Assume FF verifies 1) 2) 3). Let I' = reflexion(I,W). If IF 
does not satisfy 4), the set F(cut) = <^ a ei(W\cut)C(o-) does not cut W in 
an open set. Thus is contained in one side of the hyperplane IF. The set 
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I(W\cut) is left stable under the procedure reflexion. Clearly, the other cone 
F' (facet) = r\ a £r(w\facet)C(cr) is on the other side of the hyperplane W. 
Thus the set I 1 is not feasible. 

The following result justifies the difficulty of finding the combinatorial 
chamber that contains an input vector: 

Proposition 37 Let A be an integral matrix. Let a vector b in the cone 
C(A) generated by the columns of A and a list F of simplicial cones with 
rays in the columns of A such that all elements of F contain b. Deciding 
whether F includes all simplices that contain b, i.e. whether F determines 
the combinatorial chamber that contains b, is NP-hard. 

Proof: One well-known NP-complete problem is that of given a complete 
graph with positive integral weights on the edges to decide whether there is 
a hamiltonian tour of cost less than j3. We will explain now why this can 
be transformed of the problem of deciding whether a list simplicial cones is 
already enough to determine a chamber. 

We will use a theorem by K. Murty (see Theorem 2.1 in ^01 ) : Consider 
a complete bipartite graph K n ^ n . Orient the edges all in the same direction 
and assign excess 1 to the tail nodes and -1 to the head nodes of each arc. It is 
well known that the associated Network polytope is the famous Birkhoff-Von 
Neumann polytope of doubly stochastic matrices we saw in the introduction. 

2 

This polytope is embedded in R n and the coordinates are in correspondence 
with the arcs of the bipartite network. The associated network matrix has 
rank 2n — 1, 2n rows and n 2 columns one per arc in the network and we 
label them (1, 1), (1, 2), . . . , (n — 1, n), (n, n). 

Extend the above network matrix by adding a row of costs, where Ci», 
i ^ j, is the cost to go from i to j, except for the entry associated to the arc 
i, i where one can put a huge integer value M, much larger than the sum of 
the n largest Qj's. On the righthandside of the matrix equation we add an 
entry of value (3. Written in terms of equations we have 
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%ij > 0, for alii, j. 

This system has now rank In. The important point is: If the set of 
columns {(1, 1), (2,2), . . . , (n,n), (ixJx), {i 2 , h), (i n ,3n)} defines a sim- 
plicial cone containing the vector b = (1, 1, 1, ... , —1, —1, —1,(3), then (ii, ji), . . . , (i n ,jn) 
must be a traveling salesman tour with cost less or equal to (3. Thus if we 
take as F the set of all simplicial cones of bases that do not use all columns 
{(1, 1), (2, 2), . . . , (n, n)} and contain b = (1, 1, ... , —1, —1, /3), the remaining 
job of deciding whether any other cone contains the vector b is then at least 
as hard as the solution of the traveling salesman problem. 

To conclude this section it is worth mentioning that one can abstractly 
apply reflexions to the non-essential walls satisfying 1) 2) and 3). The 
interior of the resulting "chamber" may actually have empty interior in that 
case and thus is not useful for us here. Nevertheless this phenomenon plays 
an important role in the theory under the name of virtual chambers. In fact, 
there is another characterization of the chambers using the triangulations 
of the Gale diagram of the original vectors (see [23] for an introduction to 
Gale diagrams and triangulations). 

Lemma 38 (See [5, 10J) The face lattice of the chamber complex of a vec- 
tor configuration A is anti-isomorphic to the face lattice of the secondary 
polyhedron of the Gale transform of A, A. The vertices of the polyhedron 
are the regular triangulations of A. 

Thus generating the chambers of a network cone is the same as gener- 
ating the distinct regular triangulations of the Gale diagram of an extended 
network matrix. Such calculations can be also be done using the software 
topcom. 

4 Computational Experiments 

Now we present some computational experiments. All experiments were 
done in a 1 GHZ pentium computer running Linux using Maple 7. All 
our software is available at www.math.ucdavis.edu/~totalresidue We 
present our experiments in three tables. We begin with Table Q and Ta- 
bic |3 that deal with Kostant's partition function, this is the case of acyclic 
complete graphs. As we saw in Lemma ^ all other networks can be em- 
bedded into this case. We did examples in the cases of K4 (A3 ), (A^) 
in the first table and in the second table we have bigger examples for the 
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cases Aq Ay , Ag, Ag and Af . We show computation times in both tables 
and Table |21 also shows the cardinality of the special permutation sets. The 
computations show that the total residue method is faster than brute force 
enumeration and the current implementation of software LattE by one 
or two orders of magnitud. LattE, on the other hand, is the only software 
that deals with arbitrary rational convex polyhedra. 

As it is clear on the two first tables, the computation time does not 
increase significantly when the weights on nodes are very large. In contrast, 
computation time becomes quickly very large, when the number of nodes 
on the graph is growing. In the second table it is evident that for a fixed 
number of nodes, time of computation depends strongly of the cardinality 
of the set Sp(a), i.e. the signs of weights on the nodes (when all weights are 
positive, except the last, the cardinality of Sp(a) is 1). 

Let us stress that one of the features of our method is that it can directly 
compute the polynomial k$(a) giving the number of lattice points in the 
polytope P(&, a) in the chamber determined by a. In particular, the Ehrhart 
polynomial of the polytope P(cp,a), i.e. the function t h-> k$(ta) is also 
computed easily from our algorithm. For example, corresponding to the 
first line of Table [21 

k A + (t, 2t, 3t, At, 5t, -15t) = 1^60 ( 6 * + W + 4 ^ + 3 X* + 2 X* + *) x 

(6492U 5 + 233897t 4 + 307649t 3 + 184639t 2 + 50574t + 5040) 

which was computed in 0.55 seconds. In contrast, the polynomial function 
A;$(oi, a2, 03, 04, 05) (with 05 = —(01 + 02 + 03 + 04)) in the chamber chamber 
{ai > 0, 02 > 0, 03 > 0, 04 > 0} is computed in 0.48 seconds. 

The Ehrhart polynomials for the second, third and fourth examples 
in Table i.e. k A + (21128 * t, 45716 * t, 79394 * t, -76028 * t, -31176 * 
t, 66462 * t), k A + (82275t, 33212t, 91868t, -57457t, 47254t, -64616i, 94854£), 
and k A + (31994t,'-12275£, 55541t, 72295t, 26697i, -3212t, -38225i, 6916t), were 
computed in 1.36 seconds, 18.54 seconds, and 93.36 seconds respectively. It 
is also amusing to check the program on the value of the Kostant parti- 
tion for A^ on the vector a = [1, 2, 3, 4, • • • , r, — r(r + l)/2]. As proven by 

Zeilberger [2^, this value is Y\ r i= i nu+i)\ • 

The last table is dedicated to 4 x 4 transportation matrices. In the case 
of transportation polytopes, i.e. complete bipartite graphs. Here we also 
able to compare our speed to the special purpose C ++ program written 
by Beck and Pixton [^. Both LattE and Beck-Pixton's software are faster 
than our Maple implementation, with Beck-Pixton's significantly so, but it 
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Weights on nodes 


# of flows 


sees 


[6, 8, -5, -9] 


223 


0.1 


[9, 11, -12, -8] 


330 


0.1 


[1000, 1, -1000, -1] 


3002 


0.009 


[4383, -886, -2777, -720] 


785528058 


0.1 


[4907, 2218, -3812, -3313] 


20673947895 


0.1 


[47896, 30744, -46242, -32398] 


19470466783680 


0.01 


[69295, 62008, -28678, -102625] 


179777378508547 


0.1 


[3125352, 6257694, -926385, -8456661] 


34441480172695101274 


0.01 


[6860556, 1727289, -934435, -7653410] 


91608082255943644656 


0.1 


[12, 8, -9, -7, -4] 


14805 


0.081 


[125, 50, -75, -33, -67] 


6950747024 


0.020 


[763, 41, -227, -89, -488] 


222850218035543 


0.019 


[11675, 88765, -25610, -64072, -10758] 


563408416219655157542748 


0.011 


[78301, 24083, -22274, -19326, -60784] 


1 108629405144880240444547243 


0.029 


[52541, 88985, -1112, -55665, -84749] 


3997121684242603301444265332 


0.010 


[71799, 80011, -86060, -39543, -26207] 


160949617742851302259767600 


0.010 


[45617, 46855, -24133, -54922, -13417] 


15711217216898158096466094 


0.21 


[54915, 97874, -64165, -86807, -1817] 


102815492358112722152328 


0.060 


[69295, 62008, -28678, -88725, -13900] 


65348330279808617817420057 


0.010 


[8959393, 2901013, -85873, -533630, -11240903] 


6817997013081449330251623043931489475270 


0.010 


[2738090,6701290, -190120, -347397, -8901863] 


277145720781272784955528774814729345461 


0.010 


[6860556, 1727289, -934435, -818368, -6835042] 


710305971948234346520365668331191134724 


0.060 



Table 1: Testing for the complete graphs K4 and K§. Time is given in seconds. Excess vectors are in the first 
column. 



Weights on nodes 


# of flows 


sees 


\Sp(a)\ 


[1,2,3,4,5, -15 ] 


5880 


0.02 


1 


[21128,45716,79394,-76028,-31176,66462,-105496] 


58733548560911702671 
16780821466940568432 
553474831987566395925 


0.22 


8 


[82275,33212, 91868, -57457,47254,-64616,94854,-227390] 


22604049468113537772 

228176193404009135 

6424181 


2.14 


26 


[31994,-12275, 55541, 72295,26697,-3212,-38225,6916,-139731] 


11446847479255704222 

87042245223206779226 

01568734727431018393 

069006356672309031382 

51984519069399479632 

6644137066000 


7.94 


24 


[12275,55541, 72295, 26697,-3212,-38225,6916,92409,9528, -234224] 


12970047729476531166 
58326881685949118367 
16319862924094634125 
27856414458487356258 
66474206451882923253 
41990044115208492747 
58896993761880000897 
382293730 


21.31 


16 


[1,2, 3,4,5,6,7,8,9,10, -55] 


38883505145515430400 


5 


1 


[46398,36794, 92409,-16156,29524,-68385,93335,50738,75167, -54015, -285809 ] 


20889867895116832060 
28578373441423712122 
50684806890637191792 
33590765780756053509 
92237184823590262176 
29560725791309259479 
21077842421668832691 
54404688022155977982 
34585056426719876125 
028873152 


2193.23 


322 



Table 2: Testing for complete graphs K n with n = 6, 7, 8, 9, 10, 11. Time is given in seconds. 



must still be emphasized that our calculations for transportation polytopes 
makes use of the fact that they are embedded inside the complete graph for 
large enough number of nodes. For example the case of 4 x 4 transportation 
polytopes is treated via the complete graph K$. The same kind of embedding 
can be done for other networks. 

If we consider the case of 4 times 5 matrices with weights on nodes 
[3046, 5173, 6116, 10928], [182, 778, 3635, 9558, 11110], the number of lattice 
points is 23196436596128897574829611531938753 calculated in 11.15 sec- 
onds. The number of special permutations for this vector is 540 while the 
number of vertices of the corresponding polytope is 912. These same exam- 
ple takes 7.8 seconds in LattE and 0.1 seconds in Beck-Pixton program. 

Ehrhart polynomial fc]> 45 ((3046 * i, 5173 * i, 6116 * t, 10928 * t, -182 * 
t, -778 * t, -3635 * t, -9558 * t, -11110 * t) is computed in 30.72 seconds. 

If we consider the case of 5 times 5 matrices with weights on nodes 
[30201, 59791, 70017, 41731, 58270], [81016, 68993, 47000, 43001, 20000] , the num- 
ber of lattice points is 

24640538268151981086397018033422264050757251133401758112509495633028, 

which we computed in 23 minutes. The number of special permutations 
needed is 9572 while the number of vertices of the corresponding polytope 
is 13150. This example took 20 minutes with LattE and just 4 seconds with 
Beck-Pixton program. 

Transportation polytopes were treated by Beck and Pixton [I] in a special 
purpose C ++ program dedicated for this particular family of flow polytopes. 
Their computation is also via residues and is the fastest at the moment. It 
is important to remark that their use of residues is quite different from ours; 
our main theorem can be thought of as a multidimensional analogue of the 
fact that sums of the residues of a rational function on Pi(C) is zero. It 
is to be expected that in a forthcoming C ++ implementation the timings 
discussed here will be considerable faster than those from this preliminary 
Maple implementation. Besides obvious implementation speed ups, the ideas 
presented in this paper could still be improved when the total residue method 
is applied directly to the bipartite graph, not as a subnetwork of K n . 
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